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Abstract — We show that given n and fc, for q sufficiently large, 
there always exists an [n, k\ q MDS code that has a generator 
matrix G satisfying the following two conditions: 
(CI) Sparsest: each row of G has Hamming weight n — k + 1; 
(C2) Balanced: Hamming weights of the columns of G differ 
from each other by at most one. 

I. Introduction 

We study the existence and provide a construction of a 
sparsest and balanced generator matrix of Maximum Distance 
Separable (MDS) codes. A generator matrix is the sparsest 
if it contains the least number of nonzero entries among all 
generator matrices of the same MDS code. A generator matrix 
is balanced if every column contains approximately the same 
number of nonzero entries. More specifically, we require that 
the number of nonzero entries in each column differs from 
each other by at most one. 

Apart from being of theoretical interest, our study on 
balanced sparsest generator matrices for MDS codes was 
motivated by its application in error correction for sensor 
networks. Suppose n sensors, Si, .. ., S n , collectively measure 
k conditions x\,. . . ,Xk, such as temperature, pressure, light 
intensity, etc. Let x = (xi, . . . ,Xk), where Xi £ ¥ q for each 
i = 1, . . . , k (¥ q is a finite field of q elements). These sensors 
transmit the information they collected to a base station, which 
is a data collector. Furthermore, each sensor performs some 
encoding on the information it has, before transmitting the 
information back to the base station in the following way. Let 
G be an k x n generator matrix of an [n, k, d] q error-correcting 
code. Sensor Si transmits the scalar product of x and column 
i of G to the base station. It is well known in classical 
coding theory that this coding scheme allows the base station 
to retrieve x when at most [^"'"J sensors transmit wrong 
information. Moreover, the base station can also identify 
the malfunctioned sensors. For each sensor Si, only those 
conditions corresponding to nonzero entries of column i of G 
are involved into encoding. So it is sufficient for Si to measure 
only such conditions. Thus, if G is sparse then in average, 
each sensor only needs to measure a few among k conditions 
in order to achieve the desired error correction capability. On 
top of that, if columns of G have approximately the same 
number of nonzero entries then the sensors are required to 
measure approximately the same number of conditions. This 
balance guarantees an even distribution of workload among 
sensors, which is an important criterion for sensor networks 
where energy saving is a critical issue. 



In fact, any error-correcting code can be used in the 
aforementioned scheme for sensor networks. We choose to 
study MDS codes first because their structure, especially their 
weight distribution, is well studied (see, for instance fl] Ch. 
11]). Moreover, they have optimal error-correcting capability, 
given the length and the dimension. We prove that over a 
sufficiently large field, there always exists an MDS code that 
has a balanced and sparsest generator matrix, which is ideally 
suitable for the above encoding scheme for sensor networks. 

Necessary notations and definitions are provided in Sec- 
tion __] We state and prove our main result in Section ___\ 

II. Preliminaries 

We denote by ¥ q the finite field with q elements. Let 
[n] denote the set {1, 2, . . . , n}. The support of a vector 
u = (u\, . . . ,u n ) £ ¥™ is defined by supp(u) = {i £ [n] : 
Ui ^ 0}. The (Hamming) weight of u is |supp(w)|. We can 
also define weight and support of a row or a column of a 
matrix over some finite field, by regarding them as vectors 
over that field. Apart from Hamming weight, we also use 
other standard notions from coding theory such as minimum 
distance, linear [n, k] q and [n, k,d] q codes, MDS codes, and 
generator matrices (for instance, see flT)). 

For a matrix G = (gij) £ ¥ q xn , the support matrix of G, 
denoted supp(G), is a k x n binary matrix TVf = (rriij) where 
rriij = if gi.j = and m^.j = 1 if gij ^ 0. Let M = 
(rriij) be a k x n binary matrix. We denote by var(Af) = 
(vi.j) the matrix obtained from M by replacing every nonzero 
entry m^j = 1 by ^j, where £ij's are indeterminates. More 



formally, Vi 



if TYli 



and Vi 



if in,,, 



1. 



We also denote by gr(M) the bipartite graph Sf = (f, S) 
defined as follows. The vertex set ~y can be partitioned into 
two parts, namely, the left part L = {£%, ...,£).}, and the right 
part R = {ri, . . . , r„}. The edge set is 



*={(ti,Tj) : i £ [k], j £ 



For any k x n matrix TV, we define f(N) = J| p dct(P), 
where the product is taken over all (?) submatrices P of order 
k of TV. 

III. Main Result 

A sparsest generator matrix of an [n, k] q MDS code would 
have precisely n—k+1 nonzero entries in every row. Moreover, 
if it is balanced, then each column contains either I fc ( fl ~ fc+1 ) | 
or [ n ~ + ' 1 nonzero entries. Hereafter, we often use Ri, 



i G [k], and Cj, j G [n], to denote the supports of row i and 
column j, respectively, of a k x n binary matrix iVf. Note that 

i?i C [n] and Q C [k]. 

Lemma 1. Let M = (m^) be a kxn binary matrix. Suppose 
that each row of M has weight n — k + 1. Then M is the 
support matrix of a generator matrix of some [n,k] q MDS 
code over a sufficiently large field ¥ q (q > (?Zi) ) if and only 
jf/(var(M))#0. ' ' 

Proof: Suppose M = supp(G), where G = (gij) 
is a generator matrix of some [n, k\ q MDS code. Due to 
a well-known property of MDS codes (see [1, p. 319]), 
every submatrix of order k of Q has nonzero determinant. 
Therefore, /(G) ^ 0. Note that /(var(Af)) can be regarded 
as a multivariable polynomial in F g [. . . , ■• .]. Moreover, 
since M = supp(G), we deduce that /(G) can be obtained 
from /(var(Af)) by substituting £y by g^j for all i,j where 
gij ^ 0. As /(G) ^ 0, we conclude that /(var(M)) ^ 0. 

Now suppose that /(var(Af)) ^ 0. Note that each column 
of var(Af ) belongs to precisely (2Z-\) submatrices of order k 
of var(M). Hence the exponent of each in /(var(TVf)) 
is at most Since /(var(M)) ^ 0, by (21 Lemma 

4], if g > (fcZi) men th ere exist gi j e F 9 (for i,j 
where m^j = 1) so that /(var(iW))(. . . ,gij, ■ ■ ■ ) ^ 0. Let 
G = (gij) (for i,j where mi j =0 we set gi j = 0). Since 
/(G) ='/(var(M))(..., 5iJ ,.'..) ^ 0, again by ffi p. 319], 
we deduce that G is a generator matrix of an [n, k] q MDS 
code. Therefore, each row of G has weight at least n — k + 1, 
due to the Singleton Bound (see HI p. 33]). Since each row 
of M also has weight n—k + 1, we deduce that gij ^ 
whenever rriij = 1. Therefore, M = supp(G). ■ 

Lemma 2. Let M = (m^j) be a k x n binary matrix. 
Then /(var(JVf)) ^ if and only if every bipartite subgraph 
induced by the k left-vertices and some k right-vertices in 
gr(M) has a perfect matching. 

Proof: Let = gr(iW). Each submatrix P of order k 
of var(iVf) corresponds to a bipartite subgraph J#p induced 
by the k left-vertices and some k right-vertices in . In the 
literature, P is usually referred to as the Edmonds matrix of 
Jfp. It is well known (see J3] p. 167]) that a bipartite graph has 
a perfect matching if and only if the deteminant of its Edmonds 
matrix is not identically zero. Hence the proof follows. ■ 

Lemma 3. Let M = (m^j) be a k x n binary matrix. Then 
every bipartite subgraph induced by the k left-vertices and 
some k right-vertices in gr(M) has a perfect matching if and 
only if 

| Ujgj Cj | > \ J\, for every subset J C [n], \J\ < k. (1) 
Proof: Let S? = gr(Af). Each submatrix P of order k of 
var(TW) corresponds to a bipartite subgraph fflp induced by 
the k left-vertices and some k right-vertices in <S . The lemma 
follows by applying Hall's marriage theorem to each of such 
subgraphs of gr(M). ■ 

Lemma 4. Let M = (rriij) be a k x n binary matrix. The 
condition (fJ]) is equivalent to 



| U e j Ri\ >n-k+ |/|, for every subset ^ / C [k]. (2) 
Proof: Suppose that (Q]) holds and that there exists a 
nonempty set I C [k] satisfying 

\u ieI Ri\<n-k+\I\-l. (3) 

We aim to obtain a contradiction. The condition (O is equiv- 
alent to 

|n ie/ 7^| >fc-|/| + i. (4) 

Hence there exists a set J of k — \I\ + 1 columns of M that 
satisfies | flj e j Cj\ > \I\. Equivalently we have 

\u j ejC j \<k-\I\<k-\I\ + l = \J\. (5) 

We obtain a contradiction between ([T]i and (0. The "only if 
direction can be proved in a similar manner. ■ 

Lemma 5. Let M = (mij) be a kxn binary matrix. Suppose 
that each row of M has weight n — k + 1. Then M is the 
support matrix of a generator matrix of some [n, k] q MDS 
code over a sufficiently large field ¥ q (q > (tZi) ) if and only 
if (0 holds. 

Proof: The proof follows from Lemma [T]|4] ■ 

We present below our main result. 

Theorem 6 (Main Theorem). Suppose 1 < k < n and q > 
(feZi)- Then there always exists an [n, k] q MDS code that has 
a generator matrix G satisfying the following two conditions. 
(CI) Sparsest: each row of G has weight n — k + 1. 
(C2) Balanced: column weights of G differ from each other 
by at most one. 

By Lemma |5j to prove Theorem |6] we need to show that 
there always exists a k x n binary matrix M satisfying the 
following properties 

(PI) each row of M has weight n — k + 1, 
(P2) column weights of M differ from each other by at most 
one, 

(P3) | U ie i Ri\ >n-k+\I\, for every subset [k], 
where Ri denotes the support of row i of M. 

We prove the existence of such a binary matrix by designing 
an algorithm (Algorithm 1) that starts from an initial binary 
matrix which satisfies (PI) and (P3). In each iteration, the 
matrix at hand is slightly modified so that it still satisfies (PI) 
and (P3) and its column weights become more balanced. When 
the algorithm terminates, it produces a matrix that satisfies 
(PI), (P2), and (P3). 

Observe that it is fairly easy to construct a binary matrix 
that satisfies (PI) and (P2), using the Gale-Ryser Theorem (see 
Manfred [4]). However, (PI) and (P2) do not automatically 
guarantee (P3). Indeed, the matrix P given below satisfies 
both (PI) and (P2). However, (P3) is violated if we choose 
I ={1,2,3}. 
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Let M be any k x n binary matrix that satisfies both 
(PI) and (P3). For instance, we can shift the vector 
(1 1 ■ ■ ■ 1 • • • 0) k times cyclically to produce k rows 



n — k+1 

of such a matrix as below. 



M 



A 







o 





\0 ••• 1 1 1 1 ••• 1/ 
The Algorithm 1 takes M as an input parameter. 



Algorithm 1 

Input: n, k, M; 
Initialization: M := M; 
i: repeat 

2: Let max and min be the maximum and minimum 

weights of columns of M ; 
3: if max — min < 1 then 
4: Return M; 
5: end if 

6: Find two columns j max and j m i n that have weights 

max and min, respectively; 
7: Find a row i s satisfying TOi sJmax = 1 and 

TO j s ,imin — and moreover, if we set m.i sJmax := 

and J7ij s ,j min := 1 then M still satisfies (PI) and 

(P3); 

8: Swapping: set TOi S) j max = and mi t j min := 1; 
9: until max — min < 1; 



Due to space constraint, we have prepared a separate note 
at 1 5 1 with an example to demonstrate the algorithm. 

Lemma 7. Suppose in every iteration, Algorithm 1 can always 
find a legitimate row described in Step 7. Then the algorithm 
terminates after finitely many iterations and returns a matrix 
satisfying (PI), (P2), and (P3). 

Proof: At a certain iteration, let A = max — min. After 
swapping the two entries m.j sJmax and Trii et j min , the weight 
of column j max is decreased by one whereas the weight of 
column j m i n is increased by one. Therefore, after at most 
[n/2\ iterations, A is decreased by at least one. Hence, the 
algorithm must terminate after finitely many iterations. The 
ouput matrix obviously satisfies (PI), (P2), and (P3). ■ 

Lemma 8. In every iteration of Algorithm 1, a row i s as 
described in Step 7 of the algorithm can always be found. 

Since column j max has a larger weight than column j mm , 
there always exists at least one row i s where rn^j^ 



1 



and mi 



0. Obviously, swapping m. 



and m<i 



does not make M violate (PI). The stricter criterion is that 
M must still satisfy (P3) after the swap. We need a few more 
auxiliary results before we can prove Lemma [8] 

Suppose at a certain iteration, we choose some columns 
Jmax and j m i n that have maximum and minimum weights, 



respectively. Without loss of generality, we assume that the 
first t rows are all the rows of M satisfying the property that 
each of them has a one at column j max and a zero at column 
jmin- In other words, assume that 

{i G [k] : nii j max = 1 and m, . Jmin = 0} = [<]. 

Since max — min > 2, we have t > 2. 

Suppose, for contradiction, that none of these t rows satisfy 
the condition in Step 7 of Algorithm 1. Let M (i) , i G [t], be 
the matrix obtained from M after swapping the two entries 
TO ijmax and m,i Jmin . Then i G [t], does not satisfy (P3). 

Since M satisfies (P3) and the only difference between Af ^ 
and M is the row i, the set of rows of Af^ that violates the 
condition (P3) must contain row i. Therefore, for each i € [t], 
there exists a set U C [k], i ^ 1%, such that {i} U U is a set 
of rows that violates (P3) in For our purpose, for each 

i£ [t], we choose U to be of minimum size among those sets 
that satisfied the aforementioned requirement. Since for each 
i G [t], l-R^I = \Ri\ = n - k+ 1, we deduce that U ^ 0. 

Let Rr denote the support of row r of M^ l \ i G [t], r G 
[k]. Note that R r denotes the support of row r of M , r G [k]. 
For simplicity, we use Rj to denote the union U re _r-Rr for 
any subset / C [k]. Since {i} U U is the set of rows of 
that violates (P3), for every i G [t] we have 

\R{luu \<n~k+ \{i} U h\ - 1 = n - k + \h\. (6) 
Lemma 9. For all G [t], the following statements hold 
,(<) = \Ri>, ifi' ^ i, 

\(Ri> \ {jmax}) U {jmin}, if %' = h 
C) jmax t Rl % , d) 

= n-k + \L\. 



a) Rl 



b) =R h , 



e) iil v , f) |^. }UJi , 



Proof: Proof of a). Note that all the rows of M^> except 

(i) 

for the row i are the same as that of M. Therefore, RY,' = R4, 
if i' =^ i. As row i of is obtained from row i of M by 



swapping vtii 



1 and m, 



0, we deduce that 



Rf = (Ri \ {jmax}) U {jmin}- 



Proof of b). By definition of i £ Therefore, using Part 
a), we conclude that i?^ = J?j 4 . 

Proof of c). Suppose, for contradiction, that j max G J?j 4 . Due 
to Part a) and b), we have 



R 



(i) 

{i}UIi 



Rf 1 u i?[f 



= {{Ri \ {jmax}) U {jmin}) U R U 
= ((Ri\ {jmax}) URl,) U {j mi „} 
= (i?. t U i?,/ s ) U {jmin} 2 ^{i}UZf 

As iVf satisfies (P3), we have 



\R 



(i) 

{i}Uh 



> |%}u/J >n-*r+|{i}U/i| = n-*r+|Ji| + l. 



This inequality contradicts ©. 

Proof of d). Suppose, for contradiction, that j„ 

by Part a) and b) we have 



R h . Then 



R 



{i}U/ ( = ((^ \ {Jmax}) U {jmin}) U D {j min } U i2 7< 



3 



Therefore, using the fact that M satisfies (P3), we deduce that where the last equality can be explained by the fact that 

\R% UIi I > 10'min} UR h \ = l+ \R U \>l+n-k+\L t \ 
This inequality contradicts (0. 

Proof of e). Note that j max G i?;. However, by Part c), j max 
Rj .,. Hence, i ^ly . 
Proof of f)- Using Part a) we have 



R{i}uk Q R {i}ui t \{ty Indeed ' we have 



\R 



(i) 

{i}Uh 



B^UR?)] = I^UiZ/J > \R h \>n~k+\h\, 

(7) 

where the last inequality comes from the fact that M satisfies 
(P3). Combining © and (O, the proof of f) follows. ■ 

Lemma 10. For all i, i' G [t], i ^ £', it holds that IiCiIi> = 0. 

Proof: Without loss of generality, we prove that I\ n/2 = 
0. Suppose, for contradiction, that there exists I G 1% f)l2- We 
first present three claims, which are used later in this proof. 
Claim 1 : For £ = 1 , 2 we have 



\R 



(i) I 
{i}Uh\W l 



n-k + \Ii\ 



and 



Rf 



/?(') c R {i) 



(8) 



(9) 



Proof of Claim 1: Indeed, because of the minimality 
of Jj, the set {i} U (I* \ t) does not violate (P3) in M (i) . 
Therefore, 



II? 



{i}UJi\W 



I > n - k+ \{i} U (Ii \£)\=n-k + \Ii\. 



On the other hand, R 



G J? 



which also has 



cardinality n — k + |7j|, due to Lemma |9]f). Therefore 



i? 



»U/AW 



i? 



(i) 

{i}uJi 



and 



\R 



W 1 

{i}UJi\Wl 



|i? 



(i) 

{i}U/i 



We also deduce that i?< l) C 



-k + \Ii\. 

By Lemma |9] a), we 



1 ^ "{i}UJ 4 \{i} 

have i?^ = i?^ . Thus we complete the proof of Claim 1. ■ 
Claim 2: Let K = (h\ {£}) n (I 2 \ {£}). Then for £ = 1, 2, 
the following holds 

|-R{i}UJ 4 \W \R{t}UK \ {jmax}| < l-K/nu/A^} I - I- R {^}ux|- 

(10) 

Proof of Claim 2: Using Lemma [9] a) and b), we have 

R{i}Uli\{£} = Ri U ^MOT 



= i?, U i? 



(<) 



((firu{w)\{i m m})u4; ) vw 



C{j ma 4U(i?f Ui^ w ) 

= {imax} U -Rj.jj.u/.^f}- 

Therefore, i? W u/ATO \ Unax) ^ 4}u/,\W Hence 

\(R{i}UIi\{£} \ R{e}UK) \ {jmax}| = |(-K{»}UI 4 \{«} \ {jmax}) 

\ R{£}uk\ 



< |i? 
= li? 



UK I 



(0 



-li? 



R t ^ %}u/Am 

due to ||9). Moreover, X C ij\{£}. Hence the aforementioned 
inclusion holds. We complete the proof of Claim 2. ■ 
Claim 3: If h n I2 = {^} then for £ = 1, 2, we have 

|i?MU/Am \i?A0max}| < " 1. (11) 

Proof of Claim 3: Applying (TTOb with if = 0, we obtain 

\R{t}Ul,\W \ R e\ O'max}! < |-R{i}uJA{<}l ~ 

i (n-k+\Ii\)-(n-k + l) 
= \Ii\-l. 

We complete the proof of Claim 3. ■ 
The remaining of the proof of Lemma [10] is divided into 
two cases. Our goal is to obtain contradictions in both cases. 
Case 1: h n h = {£}■ 
We aim to show that 



\R{i,2}uhui 2 \<n-k+ |{1, 2} U/iUJ. 



2 ■ 



(12) 



This is a contradiction of our assumption that M satisfies (P3). 
Firstly, since I\ H I2 = {£}, we have 



I'll 



1. (13) 



n - k + \{l, 2} U h U h\ = n - 
Secondly, we consider 

i?{i,2}u/iu/ 2 — R e u (i?{i}u/AOT \ Re ) u i R {2}ui 2 \{i} \ Re) 

= Re U {jmax} U ((-R{1}U/AW \ ^) \ {imax}) 
U ((i?{2}u/ 2 \{^} \ i?£) \ {jmax}). 

Therefore, 

|-R{l,2}U/iU7 2 I ^ \ R i\ + 1 + |(^{l}U/i\W \ i?f) \ {jmax} 
+ l(i?{2}U/ 2 \{^} \ Rt) \ {jmax} 

G3 

< (n-k + l) + l + (\h\-l) + (\I 2 \-l) 
= n-k+ \h\ + \h\- 

(14) 

Combining ( fT3l and ( fl~4T >, we obtain ( fT2l . We complete the 
analysis of Case 1. 

Case 2: (h \ {£}) n (i 2 \ {^) - ^ 7^ 0- 
We aim to prove that 



and 



Re C i?i<- , 



R K \=n-k+ \K\ 



(15) 
(16) 



If both © and (O hold then 



li? 



{i}uk 



= \R K \ = n-k+\K\<n-k+ \{£} U K\ 



which contradicts our assumption that M satisfies (P3). 
Let 8 = \R t \ R K \ > 0. As M satisfies (P3), let 



\Rk\ 



\K\+E, 



4 



where e > 0. Then 

\R{H}uk\ = \Rk\ + \Ri\R K \ = n-k+\K\ + e + 5. (17) 
We have 

-R{i,2}u/iU/ 2 = R{£}uk U (R{i}uh\{e} \ R{£}uk) 
u {R{2}ui 2 \{t} \ R{£}uk) 

= R{£}UK U {jmax} 

U ({R{l}UIi\{£} \ R{£}Uk) \ {jmax}) 
U ((R{2}UI 2 W} \ R{£}Uk) \ {jmax})- 

Therefore 

|-R{l,2}U7iUI 2 I < \R{£}UK U {jmax}| 

+ \{R{l}Uh\{£} \ R{£}UK) \ {jmax}| 
+ \(R{2}Ul 2 \{e} \ R{t}UK) \ {jmax}| 



GZ3 

< \R 



{i}uK\ 



1 



\R {1) I 
l n {i}u/i\{«}l 



\R 



\R 
(i) 



(2) | 
{2}U/ 2 \WI 



{e}uK I 

{i}UJf| 



(18) 



- \R{£}uk \ 4 

©G3 



IB 



(2) | 
{2}U/ 2 \WI 



1 



(19) 



(n-fc+|/i|) + (n-fc+|/ 2 |) 
-(n-k + \K\+e + 8) + l 

<n-fc+|Ii| + |/ 2 |-|A'| + l. 
Moreover, as I\ fl 7 2 = {^} U Tf, we h ave 

|{1, 2} U 7 X U I a | = 2 + |7 X | + \h\ - \ {£} U JRT| 
= |/i| + |7 2 |-|i<r| + l. 
As Af satisfies (P3), from (fT8l and ( fT9] l, we conclude that 
|B {1)2}u / lU i 2 | = n-fc+|7 1 | + |7 2 | - \K\ + l. 

Therefore, all of the inequalities in (fT8l must be equalities. In 
particular, the last equality forces e = and 5 = 0. As S = 
implies that ( fT51 ) holds and e = implies that ( fT6] l holds, we 
complete the analysis of Case 2. 

In any cases, we always derive a contradiction. Therefore, 
our assumption that there exists some £ € 7l fl 7 2 is wrong. 
Hence 7i fl 7 2 = 0. It follows immediately that Ii fl ly = 
for every i, i' £ [t], i ^ i'. ■ 

We are now in position to prove Lemma [8] which in turn 
implies Theorem [6] 

Proof of Lemma® Recall that we assume that 

{i £ [k] : m ijmax = 1 and m itjtnin = 0} = [t]. (20) 

Moreover, we suppose, for contradiction, that none of these t 
rows satisfy the second condition in Step 7 of Algorithm 1 . As 
shown by Lemma|9]c), d), e), and Lemma [Tol we can associate 
to each i £ [t] a subset 7j C [k] satisfying the following 

(51) i i Ii>, for all e [t], 

(52) Jmax £ R h , for all i £ [t], 



(53) j min £ R u , for all i £ [t], 

(54) h n h = 0, for all i, i' £ [t], i ^ i'. 

Due to (S2) and (S3), for each i £ [t], there exists a row r(i) £ 
Ii that has a zero at column j max and a one at column j m i n . 
By (SI) and (S4), r(i) ^ i' for all i, i' £ [t] and r(i) ^ r(i') 
whenever i ^ i'. 



row 1 
row 2 

row t 

row r(l) 
row r(2) 

row 



t rows 



< rows 



Along the rows in the set [t] U {r(i) : i £ [t]}, the weights 
of the two columns j max and j m ; n are the same (equal to t). 
The other rows of M, because of d20l i. must contribute at 
least as much to the weight of column j m i„ as to the weight 
of column j max . Therefore, in total, the weight of column 
jmax is not larger than the weight of column j m ; n of M. This 
conclusion contradicts the fact that max > min +2. ■ 

We now discuss the complexity of Algorithm 1. In the 
initial matrix M, the difference between the maximum and 
the minimum column weights is at most k — 1. Therefore, 
according to the proof of Lemma [7] the repeat loop finishes 
after at most (k — 1)|_§J iterations. It is obvious that all 
steps in each iteration can be done in polynomial time in n 
and k, except for Step 7. It is not straightforward that the 
verification of (P3) for a given k x n matrix can be done in 
polynomial time. However, it can be shown that by considering 
a special one-source fc-sink network (of linear size in k and n) 
associated with each matrix, (P3) is equivalent to the condition 
that in this network, the minimum capacity of a cut between 
the source and any sink is at least n. On any network, this 
condition can be verified in polynomial time using the famous 
network flow algorithm (see, for instance ||6]). Therefore, 
Algorithm 1 runs in polynomial time in k and n. We omit 
the proof due to lack of space. Interested reader can find the 
proof online at 0. 
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